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We develop the Husimi map for visualizing quantum wavefunctions using coherent states as a 
measurement of the local phase space to produce a vector field related to the probability flux. 
Adapted from the Husimi projection, the Husimi map is complimentary to the usual flux map, since 
they are identical for small coherent states. By improving our understanding of the flux operator 
and offering a robust and flexible alternative, we show how the Husimi projection can provide a map 
to the classical dynamics underlying a quantum wavefunction. We demonstrate its capabilities on 
bound systems with electromagnetic fields, as well as on open systems on and off resonance. 



I. INTRODUCTION 

The probability flux, or probability current, is intro- 
duced in quantum mechanics textbooks as a determinis- 
tic operator that can be calculated, but its connection to 
experiment is often left to the reader's imagination. The 
flux operator, whose expectation over the wavefunction 
gives the traditional flux j (r, p), is defined as 
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where m is the mass of a particle in the system, and r 
and p the position and momentum. The concept of "flux 
at a point" seems paradoxical because we say something 
about momentum while also knowing position precisely. 
This raises the question: Can the flux even be measured? 

On the other hand, probability flux vanishes on sta- 
tionary states for systems with time-reversal symmetry. 
This is a shame, since strong semiclassical connections 
between trajectory flow and quantum eigenstates lie com- 
pletely hidden in the universal value of for the flux. 
Consider the example in Fig. [T] where the strong influ- 
ence of classical orbits is seen in the scarred eigenstatepQ. 
For this bound system, the flux is always zero, but when 
it is coupled to a continuum, flux becomes useful as a 
tool for examining its dynamics, even though information 
about the dynamics clearly exists before the coupling. Is 




Figure 1: A scarred eigenstate of the stadium billiard prob- 
lem is a particle in a box with the shape shown, revealing 
the strong influence of classical orbits. The traditional flux 
provides no help: it is uniformly inside the billiard. 



it possible to bridge this gap? 

By using coherent state projections, also known as 
Husimi projections[2j, we can reveal the meaning of the 
flux operator and see how to extend it to become much 
more useful. The experimental equivalent of a flux map 
has not been discussed because it is effectively impossible 
to measure - determining the flux, even at a single point, 
requires an infinite number of measurements. Instead, we 
offer an experimentally feasible extension of the flux op- 
erator based on Husimi projections which produces iden- 
tical results to the traditional flux (Eq. [T]) in the limit of 
infinitesimal coherent states. Because the Husimi projec- 
tion is able to work away from this limit and on a wider 
variety of systems, it is able to bridge the gap between 
stationary and scattering states and identify conductance 
pathways in large transport systems. 

When many Husimi projections are sampled across a 
system, they produce a Husimi map which is a powerful 
tool for interpreting the semiclassical behavior of wave- 
functions. In addition to showing the locations and di- 
rections of classical trajectories suggested by a wavefunc- 
tion, Husimi maps can also quantify how boundaries and 
external fields affect those trajectories. Previous work 
laid out the foundation for extending the flux operator 
using the Husimi projection[3j. In this paper, we present 
a complete discussion of the results summarized there 
and demonstrate Husimi maps on a wider variety of sys- 
tems with and without external fields. We then show 
how to use Husimi maps to interpret flux through vari- 
ous open devices. 



II. MEASUREMENT AND THE FLUX 
OPERATOR 

A. The Gaussian Basis 

Several discussions connecting the flux to experimen- 
tal measurement exist in the literature [4 -6J; here we offer 
an alternative view. We begin by identifying the eigen- 
states of the flux operator and giving them a physical 
interpretation. 

When discussing uncertainty, the Dirac basis implicit 
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in Eq. [T] introduces singularities which we can avoid by 
replacing the delta functions with the Gaussian basis de- 
fined as 
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where d is the number of dimensions in the system. The 
Gaussian function becomes a delta function as a — » 0. In 
the Gaussian basis, the flux operator is 

jr ,<7 = ^ (|r , o) (r , cr\ p + p |r , a) (r , cr|) . (3) 

The eigenstates, projected onto each orthogonal spatial 
dimension z, are obtained using the eigenvalue equation 



which has a solution of the form 

\K,i) = ko,cr> +api |r ,a) . 
Using the two equations 
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Further, it is useful to find the conditions on X a 
allow Eq. [8] to be written in the form of Eq. [4j 

1 



Ao-7, — 



2m 



(A 2 ) a ;A^ = 



2ma 



(4) 

(5) 

(6) 

(7) 

(8) 
that 

(9) 



Since (Pi) a = the value of a can be determined 
and from that we obtain the two eigenvalues 
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The eigenstates take the form 



(r| K,i,±) = (r\ r ,a) ± -e { • (r - r ) (r| r ,a) , (11) 
a 

where is the unit vector along spatial dimension i. 
Eq.[TT]is a linear combination of two functions: the Gaus- 
sian (Eq.[2| and its derivative. Projection of a wavefunc- 
tion onto the first term can be interpreted as measuring 
its probability amplitude at point i*o, and projection onto 
the second term as measuring its derivative along the i th 
spatial dimension at the point ro- 

Because there are two eigenstates along each spatial 
dimension, we show in Appendix [A] that determining the 



flux expectation value on a wavefunction -0(r) equates to 
calculating two dot products according to 
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Eq. [12] can be rewritten as 

if; 
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The traditional flux operator arises from the limit a — >• 
+ , at which point the two terms in Eq. 11 become the 
delta function and its derivative, while the flux values of 
the first eigenstates become 



lim X a i ± = ±00. 
7^0+ 



(14) 



In addition, there are an infinite number of other eigen- 
states with zero eigenvalues (See Appendix |A|). As a re- 
sult, a single application of the flux at a particular point 
in space j VQj i almost always results in zero, but occa- 
sionally in an extremely large positive or negative value. 
Traditionally, measurements of the flux correspond to 
the application of the operator and averaging the re- 
sults. Performing the averaging over an infinite number 
of measurements results in the expression equivalent to 
the textbook flux. 



B. Connection to Coherent States 



The prefactors before the Gaussian states in Eq. 11 
related to the Taylor expansion 
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This suggests there may be a deep connection between 
the flux eigenstates and the coherent state, defined as 
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which is a Gaussian envelope over a plane wave e* k °' r . 
Its inner product with a generic wavefunction ip (r) is 
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Observing that the phase e* k °' r ° is arbitrary, we can 
Taylor expand the exponential function in the limit of 
ko<7 <C 1 to produce 



3 



(a) (b) (c) (d) 



d/2 



(r|r ,k ,cr) 



-(r-r ) 2 /4a 2 



x(l + ik -(r-r )) (18) 
w (r|r ,cr) +zk • (r-r ) (r|r ,aXl9) 

Note that the dispersion relation for the free-particle 
continuum is a circle with radius ko = ^ 2 ™ E ; which is 
independent of the orientation of ko. The second term in 
Eq. [19] is thus proportional to the second term in Eq. [TT] 
for ko projected along the i th dimension. The similarity 
in form between Eq.[l9]and Eq. [XT] al lows us to relate the 
flux expectation value from Eqs. |12| and [13] to coherent 
state projections as 
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The traditional flux vector is constructed from the com- 
ponents in each direction. 

Unlike the Gaussian states appearing in the flux eigen- 
states in Eq. [TT] the physical meaning of coherent states 
is straightforward: they describe a semiclassical particle 
minimizing the product of position and momentum un- 
certainties. By the well-known uncertainty relation 
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taking a —> results in coherent state measurements with 
infinite uncertainty in /c-space, and zero uncertainty in 
real space. This is the limit where the traditional flux 
operates. 



C. Definition of the Husimi Projection 

The properties of coherent states make them a suitable 
basis for expanding the flux operator to a measurable 
definition, which we call the Husimi function[2j. It is 
defined as a measurement of a wavefunction ip(r) by a 
coherent state, or "test wavepacket", written as 



Hu(r ,k ,cr;^(r)) = |(^| r , k , a)\ z 



(22) 



Weighting each of these measurements by the wavevector 
produces a Husimi vector; plotting all Husimi vectors at 
a point produces the full Husimi projection. These pro- 
jections are the sunbursts in Fig. [2j which shows Husimi 
projections for the wavefunctions 
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where ki points towards the upper-right. We show the 
wavevectors that generate each state in the white arrow 
overlay. 




Figure 2: Husimi vectors for 32 equally-space points in k- 
space at left for the two wavefunctions at right: the com- 
plex plane wave (\£a) and the cosine wave (\£b) defined in 
Eq. |23| The uncertainty for each projection corresponds to 
Ak/k = 2%(a),10%(b),50%(c),250%(d), corresponding to 
smaller wavepacket spreads (middle) and less distinction be- 
tween independent measurements (top and bottom). Above, 
we represent the coherent wavepacket spread using schematic 
circles; in general, we indicate the spread using double- arrows. 



Both wavefunctions are pure momentum states which 
are not spatially localized, and constitute the building 
blocks for the wavefunctions addressed in this paper. The 
plane wave \£a is relevant to magnetic field states dis- 



cussed in Section MB The cosine wave \I/b corresponds 
to time-reversal symmetric wavefunctions discussed in 
Sections |III A| and |III C[ Both \£a and \I/b are important 



for scattering wavefunctions presented in Section |IV A| 
which exhibit a mixture of both properties. 

Because of the large momentum uncertainty for small 
<j, coherent state projections merely reproduce the prob- 
ability amplitude |^( r )| 2 m a U directions of ko, as seen in 
Fig. [2]i. The flux emerges as a small residual which can 
be retrieved by summing each coherent state projection 
weighted by ko- We call this quantity the vector- valued 
Husimi flux, 
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/ ko|<^| 
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In Appendix [Bj we show that as a — > 0, the contribut- 
ing points in the integral over /c-space reduce to just the 
orthogonal directions. In this limit, we can write the 
Husimi flux as 



lim Hu(r ,cr;^(r)) oc V ej|(^| r , k^, a)\ 2 
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-|(^|r ,-^oe,,a)| 2 ], (25) 
where is the unit vector along the i th orthogonal di- 



rection, and we sum over d dimensions. By Eq. [20j both 
sides of Eq. [25] are proportional to the traditional flux 
measured at point tq so that 
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For larger a, reduced momentum uncertainty allows for 
substantial variation in the coherent state projections be- 
tween different directions of ko. This can be seen in Fig. [2] 
as uncertainty is reduced and uniform sunbursts (d) con- 
tract into lobes (c), and finally to unambiguous vectors 
(a-b). At all uncertainties, the absence of flux in time- 
reversal symmetric states can be interpreted as the per- 
fect cancelation of coherent state projections along each 
direction in /c-space. The equal participation of counter- 
propagating flux while absent in \I/a is evident in \I/b as 
a reflected sunburst. 

The reduced momentum uncertainty for larger coher- 
ent states also reduces spatial resolution. In the interme- 
diate regime, we can use Husimi projections to map the 
local phase space of a wavefunction. By taking snapshots 
of the local phase space at many points across a system 
for larger a, we can produce a map of the classical tra- 
jectories that correspond to a given wavefunction. These 
visualizations are known as "Husimi maps" [2j [7H9] . Like 
the traditional flux map, Husimi maps can be integrated 
over lines and surfaces to reveal the total probability flux 
current. 

To produce the Husimi map, we sample Husimi pro- 
jections along a grid in spatial coordinates, since it is 
easier to plot, straightforward to interpret, and allows for 
computing spatial derivatives (see Section |IIE ). However, 
other schemes may be preferred. In Fig. [5J for example, 
we sample along classical trajectories to emphasize the 
quantum-classical correspondence. While this paper ad- 
dresses two-dimensional systems, Husimi projections are 
equally applicable for higher-dimensional systems. 

Husimi maps also have implications for experiments 
since they could be measured in a fashion similar to 
angle-resolved photo-emission spectroscopy (ARPES), 
which is currently used to measure the dispersion re- 
lation and Fermi surfaces (for a review, see [10J). In 
the ARPES setup, a focused photon beam on a sample 
kicks off electrons in the valence band. The energy of the 
photo-emitted electrons incorporates both their bonding 
energies, which can be averaged over, and their kinetic 
energy, which depends on the angle of the beam with 
respect to the sample surface. 

The ARPES response function behaves similarly to co- 
herent state projections with ko proportional to the beam 
angle. By rotating the beam angle around the same point 
of intersection, the response in different directions pro- 
vides the momentum distribution of the wavefunction at 
that point. Perturbations from the known dispersion re- 
lation can then be inserted into Eq. [24] to obtain the flux 
expectation value. 

While a narrow beam would make it possible to mea- 
sure the flux vector at the intersection point, it will be 
difficult to distinguish the occasional large perturbation 
measurements from noise. However, wider beams would 




(d) . (e) 



Figure 3: Husimi vectors for 32 equally-space poi nts in k- 
space (a) for the double cosine waves (^c) from Eq. [27] The 
uncertainty for each projection corresponds to Ak/k = 30%. 
As the multi-modal algorithm (Algorithm [I]) loops through 
each iteration (b-e), a trajectory is matched and then sub- 
tracted from the full Husimi projection, until all major tra- 
jectories are approximated by their appropriate values. 



capture additional terms from the Taylor expansion of 
the coherent state in Eq. [19] producing more reliable 
perturbation measurements. Applying the technique at 
many points across the sample would then provide the 
Husimi map and an approximation to the flux map. 

A question arises regarding the handling of boundaries 
in the system, beyond which the wavefunction goes to 
zero. Our definition reduces the magnitude of Husimi 
projections within distance a of the boundary. When 
a coherent state interacts with a boundary, the bound- 
ary can be replaced by an image wavepacket moving 
in the opposite direction. In this case, reflections off 
the boundary amount to scattering between wavepackets 
with different wavevectors. Thus, the reduction in the 
Husimi projections near the boundaries is the result of 
wavepacket scattering, making it possible to use Husimi 
maps to compute scattering metrics along the boundary, 
such as angular deflection presented in Section [HE] 



D. Multi-Modal Analysis 

The Husimi projections in Fig. [2] reveal that even a 
single plane wave produces a range of Husimi vectors be- 
cause of the finite spread of the wavepacket. Can distinct 
trajectories intersecting at a point be distinguished un- 
ambiguously? If the dominant plane waves at a point 
are sufficiently separated in /c-space, i.e. the momen- 
tum uncertainty of the coherent state can resolve between 
them, we can retrieve their wavevectors numerically using 
Multi-Modal Analysis (MM A). This analytical tool can 
be especially useful for time-reversal symmetric systems 
where both the traditional flux and the Husimi total flux 
are identically zero. 
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Algorithm 1 Multi-Modal Analysis (MMA) 

1. A set of Husimi templates on N wavevectors {kj} is 
created for the wavef unctions ^ = e 2k * r generated by 
the M wavevectors {k^ est } . Both sets of wavevectors 
lie along the dispersion contour. Each template can be 
stored as a vector of values of length M where each 
member corresponds to the Husimi function along the 
wavevector kj. 

2. Writing the Husimi projection as the vector v, a metric 
is computed di — v • for each Husimi template. 

3. The maximum of the set {di} is determined, and both 
the wavevector k^ est and the dot product di are stored. 

4. The contribution of the trajectory with wavevector k^ est 
is determined by the re- weighted vector - d * u . 

5. The re- weighted template vector is subtracted form the 



projection, that is, v — >• v — 



6. All elements of v which are now negative are set to zero. 

7. Steps 1-6 are repeated until the metric di dips below a 
threshold. 

8. The set of vectors {c^k* est } are used to approximate 
the Husimi projection 



Figs. [3] demonstrates the MMA algorithm on the pure 
momentum state 

^c(r) = acos(ki • r) +/3cos(k 2 • r) , (27) 

where ki points towards the upper-right and k 2 points 
towards the upper- left. We set a — 1 and j3 — 0.8. In 
Fig. [3^i, the Husimi projection is shown with a sizable 
uncertainty of Ak/k = 30%. Parts b-e iterate through 
the for loop in steps 1-6 of the MMA Algorithm. At 
each iteration, the most dominant plane wave in the sun- 
burst is modeled and then subtracted from the projec- 
tion. This is repeated until all major plane waves have 
been approximated. If the dominant trajectories inter- 
secting at a point have sufficiently divergent momenta, 
not only does the algorithm do an excellent job of mod- 
eling them, it can even compute how many there are. In 
general, we stop the loop in Step 7 after a certain number 
of iterations to make clearer figures. 

On the other hand, when there are a number of trajec- 
tories of equal weight whose momenta cannot be resolved 
by the coherent state, the MMA Algorithm can produce 
unexpected results. An example of unresolved trajecto- 
ries is seen in the points sampled along the perimeter of 
Fig. [4^i and in the central regions of Figs, ^p-d and [5] In 
these cases, the MMA Algorithm approximates overlap- 
ping trajectories by first choosing their average, and then 
contributing additional trajectories on either side. 

When the traditional flux is non-trivial, as in Fig. [7] 
and magnified in Fig. [8^i, it averages over trajectories 
at each point to produce the total drift flow. For this 
reason the multi-modal analysis is able to augment the 
information provided by the flux operator since it can 
show the individual trajectories contributing to the av- 



erage. For example, in Sec. |IIIB[ we compare the drift 
flow highlighted by the flux to the classical paths and the 
multi-modal analysis for the same system. 



E. Angular Deflection 

The Husimi map makes it possible to compute other 
quantities tied to the semiclassical underpinnings of a 
quantum wavefunction. This section, and the results ex- 
amined in Section |III C[ focus on one: angular deflec- 
tion, which reveals where system boundaries and exter- 
nal fields deflect classical trajectories from straight paths 
and give rise to the shape and properties of a given wave- 
function. 

We begin by considering the Husimi function for one 
point in /c-space measured at equally-spaced points on 
a grid that covers the system. The scalar field yields a 
spatial map of the presence of an individual trajectory 
angle, and fluctuations in the map indicate points where 
classical paths deflect away from and towards the angle. 
Summing the results for all wavevectors along the con- 
tour line defined by system energy in the dispersion rela- 
tion, we can derive a measurement of angular deflection 
Qang. (r; \P) written as 

Qang. (r; *) = f D a bs. (r, k; *)kd d k. (28) 

D a bs.(r, k; \£) is the Gaussian- weighted absolute diver- 
gence of the Husimi map for wavevector k written as 
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(29) 



where we sum over the d orthogonal dimensions each as- 
sociated with unit vector e$. 

Using the Husimi map to measure angular deflection 
has close ties to its initial introduction as a measurement 
state for building phase diagrams |2J. For instance, it is 
possible to use the divergence of the Husimi map for each 
wavevector to compute the quantum analog of a state's 
Poincare map[llj. This form of the Husimi map has been 
used to examine the angle of impact against a coordinate 
along the boundarv|ll| to study chaotic behavior in sta- 
dium billiards H21 IE]. 



III. HUSIMI MAPS IN CLOSED SYSTEMS 

A. Eigenstates of the Circular System 

The circular well is an ideal system for demonstrating 
the Husimi map since their classical dynamics are simple 
and can be analytically determined. 
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Figure 4: Husimi maps (left), multi- modal analysis (middle), 
and the wavefunction (right) for eigenstates of the circular 
well (a-c) and the harmonic oscillator (d). Double-arrows at 
far right indicate the spread of the coherent state which is 
Ak/k — 10%. The states in (c) and (d) correspond to the 
classical paths in Figs. [5^t and b respectively. 




Figure 5: Quantum-classical correspondence from Husimi 
maps by sampling along classical trajectories. In part (a), 
the Husimi map for the two eigenstates in Fig. [4]>d, where 
Husimi projections are sampled along a grid. In part (b), pro- 
jections are instead sampled along classical paths that corre- 
spond to the wavefunction. Because of rotational symmetry, 
however, the wavefunction is actually created by the sum of 
many rotations of such paths, as indicated in part (c). 



The Schrodinger equation can be written in radial form 

as 

Solutions to this equation are simultaneous eigenstates 
of energy and angular momentum, and thus possess the 
good quantum numbers n (number of nodes in the radial 
direction) and m (number of angular nodes). Fig. [4^i-c 
shows three such states, the first with n = 0, the second 
with n ^> m, and the third with n « m. The Husimi 
map in each shows the clear distinction between angu- 
lar and radial components of the wavefunction, and how 
they correlate with classical paths with similar properties 
(further discussion of the classical correspondence can be 
found in Ref. [13]). 

To examine the harmonic oscillator state in Fig.[4]i, the 
Husimi projection at each point must be modifie d. For 
the circular well, the dispersion relation is hk = \j2mE, 
but due to the harmonic potential, it changes to hk(r) = 
y/2m(E-V(r)). This means that a different sweep in k- 
space must be made at each point to produce an accurate 
Husimi map. Fig. PJl shows such a state with V(r) = 
For 2 . 



The Husimi vectors in Figs. ^ align to suggest straight 
trajectories, but the vectors in Fig.[4]i do not, suggesting 
the presence of curved paths. Moreover, projections near 
the boundaries of both systems indicate that the paths 
of the circular well bounce off the boundary with a con- 
sistent and acute angle, while they graze the edge of the 
harmonic oscillator. 

In this paper, we have chosen to sample the Husimi 
projections at equally-spaced points along a grid, which 
makes it possible to compute quantities such as the an- 
gular deflection. If we instead sample along one of the 
classical paths corresponding to the state, we find a set of 
Husimi vectors which align themselves perfectly with the 
classical path. We show these two approaches in Figs. [5k 
and|5b, which correspond to the wavefunctions in Figs.^p 
and respectively. 

Each Husimi projection in Fig. ^jp contains an addi- 
tional set of Husimi vectors which do not align with 
the path. These vectors can be understood by consid- 
ering that wavefunctions for the circular well and har- 
monic oscillator actually correspond to infinitely many 
such paths rotated in space due to the circular symme- 
try of these systems, which we indicate in Fig. [5J3. The 
"cross-hatching" patterns in Fig. [5^i-b arise because two 
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Figure 6: Husimi map (left), multi- modal analysis (middle), 
and the wavefunction (right) for two eigenstates of the circular 
well with magnetic field vectors coming out of the plane. The 
magnetic field strength is set so that the cyclotron radius is 
approximately l/2(a) and l/3(b) of the the system radius. 
Double-arrows at far right indicate the spread of the coherent 
state which is Ak/k = 10%. These states correspond to the 
classical paths discussed in Fig. [7] 



rotated classical paths intersect at any point. 

Towards the center of the system, a large number of 
paths come into close proximity. Even though an in- 
finitesimal point is intersected by only two paths, the 
finite spread of the coherent state is sensitive to other 
paths nearby, giving rise to Husimi projections showing 
a large number of trajectories with similar angles. These 
points in a wavefunction can violate assumptions of the 



multi-modal analysis in Section |IID| since the different 
trajectory angles cannot be resolved by the finite spatial 
and momentum uncertainties of each Husimi projection. 
As a result, the multi-modal analysis in Figs. [4|; and|4]i 
does not produce the original paths, but their average 
and approximations on both sides of the average. 




Figure 7: The flux map, multi-modal analysis, and classical 
paths for the states represented in Fig.|6|a-b). The traditional 
flux correlates strongly with Husimi flux (Eq. |24| ) but fails to 
show the classical paths suggested by the wavefunction. Red 
circles correspond to magnified views in Fig. [8] 




Figure 8: Husimi vectors for 32 equally-space points in k- 
space are shown in grey for the double plane waves defined 
in Eq. 32 The uncertainty corresponds to Ak/k = 30%. Be- 
cause of momentum uncertainty, there is spread in the Husimi 
vectors. The flux operator (a, blue) averages over the Husimi 
vectors, while the processed Husimi map (b,red) recovers both 
underlying wavevectors. is representative of the areas cir- 
cled in red in Fig. [7] 



B. Magnetic Field 

Systems without time-reversal symmetry can also be 
studied with the Husimi technique as shown below for 
systems in the presence of a magnetic field. To prop- 
erly reflect these states, both the momentum operator in 
Eq. [T]and the momentum term zko • i*o in Eq. 24 must be 
modified to reflect the canonical transformation 



(31) 



where the magnetic potential A is defined in App. [C] 

Results for large magnetic fields, where the cyclotron 
radius is smaller than the system size (see App. JC|, are 
presented in Fig. [6] Unlike the circular well states with- 
out magnetic field in Fig. [4j the wavefunctions in Fig. [6] 
do not exhibit cross-hatching nodal patterns, but circular 
nodal patterns with complex phase arguments. Project- 
ing each phase argument onto the real axis, however, it is 



easy to see that cross-hatching nodal patterns re-emerge, 
suggesting the presence of multiple classical trajectories 
at each point. This intuition is corroborated by the cor- 
responding Husimi maps for each state, which indicate 
circular classical trajectories with radii corresponding to 
the cyclotron radius. 

In Fig. [7| the full classical paths corresponding to each 
state are depicted, and correlate strongly with the Husimi 
map with the canonical transformation. Like the circular 
well states, the presence of multiple trajectories at each 
point in Fig. [7] can be explained by the intersection of 
rotated classical trajectories that arise from rotational 
symmetry. For the state in Figs. |6^i and [T^i, we have 
artificially lifted rotational symmetry to highlight fewer 
rotated paths. 

Because the flux map averages over trajectories at each 
point, it can often fail to indicate the full classical dynam- 
ics underlying a quantum wavefunction. The left column 
of Fig. [7| which shows the flux map, integrated with a 
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Figure 9: Husimi maps for the scarred stadium billiard eigen- 
state. Each map uses a different spread of the measurement 
wavepacket. The spread is indicated by the double- arrows 
on the bottom, with relative uncertainties of Ak/k = 5% (a), 
20% (b), and 50% (c). A single Husimi projection, circled in 
red, is magnified at the bottom of each representation. 



Gaussian kernel corresponding to the coherent state used 
to generate the Husimi map, is consequently unable to 
represent the classical paths (right column), but instead 
measures and total drift flow which might be the desired 
quantity in some circumstances. In contrast, the multi- 
modal analysis in the middle column indicates the clas- 
sical paths with remarkable fidelity. 

We can appreciate the difference in detail by examining 
the areas circled in red in Fig. [7J In Fig. [8j magnified 
views from the flux operator, multi-modal analysis, and 
full Husimi map corresponding to these ares are shown. 
We model this point in the wavefunction according to the 
pure momentum state 



(r) 



Q ik 4 r 



(32) 



where k3 and k4 are indicated by the white arrows. Not 
only does the flux average the full Husimi projection, 
it also averages over the trajectories inferred by multi- 
modal analysis. As a result, the flux map in the left col- 
umn Fig.[7|can be deduced from the multi-modal analysis 
in the middle column by simply summing the vectors at 
each point. 



C. Stadium Billiard Eigenstates 

The classical dynamics of the circular stadium are in- 
tegrable while those of the Bunimovich stadium [T5j are 
chaotic. As a result, the stadium has been featured in 
many studies of "quantum chaology " fTT | H6U22] . 

Fig. [9] shows three Husimi maps for a billiard eigen- 
state. The wavelength at the energy of the eigenstate is 
much shorter than the size of the system, allowing well- 
defined scars to form, which are spawned by modestly 
unstable and infinitely rare (among all the chaotic or- 
bits) classical periodic orbits[lJ. 



For Fig. [9^i, an extended coherent state is used to gen- 
erate the Husimi map, so that many fine features of the 
wavefunction are washed out. Only the scar path (seen 
as a rotated "v" pattern in the depiction) is clearly vis- 
ible. The sharply peaked Husimi sunburst reflects both 
the low momentum uncertainty of the Gaussian used and 
the strong dominance of the periodic orbit pathway in the 
eigenfunction. 

Compare this to the Husimi map in Fig. [9^ which 
is generated by a small coherent state with larger mo- 
mentum uncertainty. Here, each Husimi projection is 
more ambiguous, and local variations in the wavefunction 
probability amplitude have a large impact on the repre- 
sentation since they are no longer smoothed over. As 
a result, the trajectories implied by the map no longer 
continue from one projection to its neighbors and ap- 
pear somewhat irregular. In general, a compromise can 
be made by choosing an intermediate momentum uncer- 
tainty, as shown in the Husimi map presented in Fig. [9)3. 
Trajectories are fairly well-resolved, and local variations 
are easy to follow. Coherent states of this size provide 
the clearest representation of semiclassical paths. 

Even at low energies, where the wavelength is compa- 
rable to the size of the system, stadium billiards pro- 
vide another perspective on the utility of the Husimi 
map. Unlike the circular system, in which the trajec- 
tories adding up at a particular point are fairly regular 
and predictable, any point in a stadium billiard eigen- 
state is rife with many unpredictable trajectories. Thus, 
the Husimi map is an ideal tool for lifting the veil on the 
underlying classical dynamics. 

In Fig. [ToJ we show Husimi maps for three eigenstates 
of the closed stadium billiard Hamiltonian. For each 
calculation, the size of the coherent state is kept con- 
stant, but because the energy of the eigenstates increases 
from top-to-bottom, the momentum uncertainty for each 
Husimi projection also increases. This is reflected in the 
clarity of the suggested classical paths at higher energy 
as well as the reduction of angular deflection in the bulk 
(which acquires small positive values in the top figure 
due to uncertainty, not because there is actual deflection 
at these points). 



To the unaided eye, the wavefunctions in Fig.[10]do not 
appear to emphasize isolated classical trajectories like the 
high-energy stadium state in Fig. [9| especially since at 
such low energies the system only accommodates a few 
wavelengths along its diameter. In the Husimi map, how- 
ever, it is quite clear that a very limited set of classical 
trajectories are largely responsible for these wavefunc- 
tions, suggesting that Husimi projections could be used 
to study the properties of low-energy scar states [Ij. 

Points with high angular deflection show which parts 
of the system boundary are responsible for the creation 
of each state, and indicate where adiabatic changes in the 
boundary conditions are most likely to affect the state [23 , 
24J. This can be imagined as a quantum force on the 
boundary. Because the size of the coherent state used to 
generate each Husimi map is kept constant, the angular 
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Figure 10: Three eigenstates of the stadium billiard system 
with Dirichlet boundary conditions at three increasing ener- 
gies, with the Husimi map (left), multi-modal analysis (mid- 
dle) and wavefunction (right). Angular deflection is indicated 
in blue, and the double arrows indicate the test wavepacket 
spread of Ak/k = 20%(a), 15%(b) and 10%(c). 



deflection penetrates into the bulk to the same extent 
for each state. However, the locations of high angular 
deflection along the boundary form a unique fingerprint 
for each state. 



IV. FLUX THROUGH OPEN SYSTEMS 

A. Sub-Threshold Resonance 

The previous section used the Husimi map to exam- 
ine the semiclassical dynamics of closed systems directly 
from their wavefunctions, providing substantial benefits 
over the usual flux operator, which vanishes for time re- 
versal symmetric systems, and averages all trajectories 
(thus missing criss-crossing trajectory paths, see Fig. [7]) 
for a magnetic field present. Moreover, the spread of the 
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Figure 11: Top: An infinite waveguide schematic with a slight 
bulge in the middle (grey). This can be modeled as two 
waveguides of different widths (blue and red). Bottom: In 
an infinite waveguide, the transmission curve has a series of 
plateaus as each transverse mode opens up (blue transmis- 
sion curve). In a wider waveguide, each mode opens up at 
lower energies (red curve). If only a small segment of the 
waveguide is widened, then sub-threhold resonances occur in 
between the energies of the narrow and wide waveguides (grey 
transmission curve) . These correlate with sub-threshold reso- 
nant states which peak in the density of states (DOS) at those 
energies (grey curve). Energy is given in units of t where 4t 
is the band edge. 



coherent state used to generate the Husimi map gives it 
the flexibility to examine dynamics at a variety of scales, 
while the flux operator is confined to the limit of in- 
finitesimal spread. In its traditional guise (Eq. 0, the 
flux operator is most often employed in scattering prob- 
lems which arise when a closed system is coupled to an 
environment. Is it possible to connect the semiclassical 
dynamics of the closed system to the open system using 
the extended Husimi flux? 

In this section, we demonstrate how the Husimi flux 
can help interpret the traditional flux and deepen our 
understanding of transport across a device. We consider 
sub-threshold resonance for a waveguide that is slightly 
widened along a short section (see inset, Fig. 11). 

In an unperturbed waveguide, transport occurs 
through transverse modes which open for transport when 
the system energy exceeds the transverse energy of the 
mode. At these energies, the transmission function ex- 
hibits distinct plateaus as seen in Fig. [TT] where the plot 
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Figure 12: The full Husimi map for the resonant state (see 
inset) is plotted with Ak/k = 100% (a) and Ak/k = 20% (b). 
The spread of the test wavepacket is indicated by double- 
arrows. A single Husimi projection (circled in red) for each 
map is magnified at right. The vector sums of each map are 
shown in Figs. |14| d. 



of the transmission for a wide (narrow) waveguide is pre- 
sented in red(blue). 

If a small section of a narrow waveguide is widened, the 
transverse energy of each mode diminishes in the wider 
section. Thus, for each mode, there is a range of energies 
bounded above by its transverse energy in the unper- 
turbed waveguide, and below by its energy in the wider 
region. In this energy range, the mode can reside in the 
wider region but cannot propagate through the narrower 
leads where it is an evanescent wave. This forces it into 
a quasi-bound state which is trapped in the wider region 
and is only weakly coupled to the environment, causing 
a striking peak in the density of states, commonly known 
as a Feshbach resonance [25J. In the quasi-bound state, 
the particle bounces vertically between the walls of the 
perturbed region and is unlikely to escape. 

At certain energies, a particle propagating in a lower 
energy mode corresponding to the narrow section inter- 
acts with the wider region and becomes trapped in the 
quasi-bound state. This causes the quasi-bound state to 
hybridize with the propagating mode and interfere with 

The 



the transmission in the device, as seen in Fig. II 



suppression of transmission appears as a pair of sharp 
dips, accounting for symmetric and antisymmetric ver- 
sions of the Feshbach resonance. Since the resulting 
wavefunction is the hybridized state which inhabits the 
system at resonance, we refer to it as the resonant state. 

We compute the wavefunction of the resonant state 
corresponding to the first transmission dip in Fig. 11 (in- 
dicated by the arrow in the transmission function) ac- 
cording to Appendix [D] Our method allows us to extract 
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Figure 13: The full Husimi map for lowest propagating mode 
wavefunction (see inset) from the waveguide in Fig. 11 is plot- 
ted in (a), at an energy well above resonance (E = 0.02745 
in arbitrary units scaled to Fig. |11| ). The uncertainty for this 
map is Ak/k = 20%. At right, a magnified view of the pro- 
jection circled in red. In (b), the Husimi flux. This is the 
mode which hybridizes with the resonance state to produce 
Figs. [12] and p] 



the pure resonant state without the second-lowest prop- 
agating mode, which is also present at these energies. 
Fig. [T2| shows the full Husimi map for this wavefunction, 
using coherent states with uncertainties of Ak/k = 100% 
(a) and 20% (b). The individual projections correspond 
strongly to the cosine- wave projections in Fig. [2] Spatial 
variations in the Husimi map decrease as the size of the 
coherent state increases, as in Fig. [9] 

The full Husimi map is indistinguishable from the 
quasi-bound state and the resonant state, which is ex- 
pected since the resonant state only slightly perturbed by 
the propagating mode. The flux of the quasi-bound state 
is zero, but exhibits characteristic vortices in the resonant 
state. Moreover, as the energy is increased across reso- 
nance, the wavefunction doesn't substantially change in 
appearance, while the flux patterns alter dramatically. 
At first these behaviors appear to contradict the Husimi 
map, but we can show that the flux patterns correlate 
with subtle changes in the Husimi maps which we can 
retrieve by adding all their vectors. 

We can begin to understand these subtle changes by 
examining the lowest propagating mode. The full Husimi 
map far away from resonance, shown in Fig.[l3]for a mod- 
erate coherent state, corresponds to the complex plane 
wave in Fig. [2] In the Husimi flux, the left-to-right flow 
appears unchanged within the central region of the sys- 
tem. The flux operator for this mode, not shown, is simi- 
lar. In contrast, the vector-sum and the flux of the bound 
state is always zero. So what happens when it interacts 
with the lowest propagating mode to produce the reso- 
nant state? 



In Fig. 14 we address this question by showing the 
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Figure 14: The traditional flux (left column) and the Husimi 
flux (right column) for the resonance state in Fig. 11 slightly 
above resonance (a, E — E ves . + 0.00005), at resonance (b, 

E — -E/res. 

) and slightly below resonance (c, E 
0.00005). The coherent state for the Husimi map corresponds 
to Ak/k = 0.2. The transmission function for this mode at 
each energy corresponds to T = 0.99(a), 0.06(b), and 0.99(c). 
Even though the full Hus imi maps at each energy are indis- 
tinguishable from Fig. 12 their vector additions (Husimi flux) 



vary substantially. Energies are in arbitrary units, scaled to 

Fig. \n\ 



across the central region which counteract the left-to- 
right flow from the leads, resulting in zero transmission 
for this mode. The second-lowest propagating mode (not 
shown), which is antisymmetric along the transverse di- 
rection, does not interact with the resonance and main- 
tains full transmission. 

At all energies, the full Husimi map shows the simple 
vertical bouncing trajectories that are identical to the 
bound state (Fig. 12), while the left-to-right flow of the 
lowest propagating mode (Fig. 13) interferes with these 
paths to produce the residual flux vortices. The classi- 
cal dynamics of the resonance therefore indicate a subtle 
shift in the overall contribution of classical trajectories 
which give rise to the resonance. Because the vertical 
trajectories can easily cancel each other out, the residual 
becomes exquisitely sensitive to the initial conditions of 
such classical paths, which are determined solely by the 
energy of the lowest-propagating mode. 

By examining this system using the Husimi flux, which 
allows us to adjust the coherent state spread arbitrar- 
ily, we can zoom into the details of the flux operator at 
small spreads and pull out to larger drift flows at larger 
spreads. Important information about the resonance can 
be retrieved at all scales, since the flow can be understood 
by the slight residuals of the full Husimi projection with 
different coherent state spreads. By adding the Husimi 
projection and the Husimi flux to the analytical toolset, 
we can examine the problem from all angles to construct 
a more nuanced and complete story. 



traditional flux, wavefunction, and the Husimi flux above 
(a), at (b), and below (c) resonance. The flux operator 
is integrated over a Gaussian kernel corresponding to a 
coherent state spread of Ak/k = 100%, and is visually 
identical to the Husimi flux with the same coherent state 
spread. 

In the flux operator, we see the characteristic vortex 
patterns which flip above and below resonance, as ex- 
pected when the bound state shifts through a phase of tt 
over resonance. Moreover, while it is clear that the pres- 
ence of the lowest propagating mode is stronger away 
from resonance, the wavefunction representation at all 
three energies are strongly influenced by the bound state. 
Similarly, probability flux is strongly localized in the cen- 
ter of the system, and it is unclear how the vortices cor- 
relate with the fact that transmission for this mode goes 
to zero on resonance. 

In the Husimi flux, however, the correlation is obvi- 
ous. Above and below resonance, vortices cancel out and 
leave behind the drift velocity of the mode. At these 
energies, the Husimi flux is quite similar to the lowest 



propagating mode in Fig. 13 , and the left-to-right flow 
extends through the semi-infinite leads, although there 
are slight changes in the central region. At resonance, 
however, the vortices no longer interfere to produce flow 
from left-to-right, but instead persist as larger vortices 



B. Transport Through Other Geometries and the 
Nature of Flux Vortices 

Because the a parameter defines the spatial spread of 
the coherent states used to generate a Husimi map, we 
can use it to reveal the behavior of the probability flux 



at arbitrary scales. In Fig. [15ft , we show the scattering 
wavefunction that acts as a mode of unit transmission in 
a large square block geometry. This geometry is changed 
from the previous subsection so that: 1) Its dimensions 
are much larger than the characteristic wavelength at the 
energies we examine, 2) the leads are shifted vertically 
from the center towards the bottom- left and upper-right, 
and 3) the center is obstructed to constrain transport 
through the central region. As a result, classical paths 
related to transport in this system must reflect off the 
boundaries many times to propagate from the left to the 
right lead. 

In Fig. [15] we show the wavefunction, a magnified view 
of the traditional flux, the full Husimi flux, and the multi- 
modal analysis for this scattering state. In the wavefunc- 
tion, nodal lines appear to fall along the 45° diagonals, 
which is corroborated by trajectories favoring those di- 
agonals in the multi- modal analysis. This arises because 
all boundary conditions are vertical or horizontal walls; 
since each mode of the unperturbed waveguide leads is 
associated with a distinct pair of trajectory angles (See 
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Figure 15: A scattering wavefunction associated with full 
transmission through the obstructed square with the wave- 
function representation (a), the Husimi flux (b) and multi- 
modal analysis (c) for a coherent state spread of Ak/k = 10%, 
indicated by the double arrows. The traditional flux from the 
part of the system indicated by the black squares is magnified 
in the inset. 



Fig. [16] ), the vertical and horizontal walls therefore re- 
flect all trajectories back onto the same pair rotated at 
45 degrees. At the energy we have selected, the pair of 
trajectory angles for the incoming mode are at perfect 
45° diagonals, so that their rotations from reflecting off 
the walls also point along the diagonals, giving rise to 
strong standing waves. 

In Fig. [TIfo , it is clear that transport occurs primar- 
ily through a narrow channel we call the conductance 
pathway with the majority of arrows pointing from the 
lower-left to the upper-right corners. By comparison, the 
full traditional flux map (not shown) is rife with vortices 
throughout the entire system, dramatically limiting our 
ability to identify overall flow. The conductance path- 
way does not have to be classical, since it is an aggregate 
phenomenon from many other classical trajectories; as 
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Figure 16: The scattering wavefunction associated with dis- 
tinct modes of a wide unperturbed waveguide at left, with the 
multi- modal analysis at right. Each mode is associated with 
a pair of trajectory angles. As the number of horizontal nodal 
lines increases, and as energy increases, these angles become 
increasingly vertical. 



a result, it is able to curve in the bulk without exter- 
nal forces, as seen in the figure. As the pathway moves 
against many other perpendicular classical paths indi- 
cated in the multi- modal analysis, pairs of vortices form 
on either side, which we show in the inset in Fig. [15] 
These vortex pairs are a direct analog to those which 
occur in sub-threshold resonance as the left-to-right con- 
ductance pathway passes through perpendicular trajec- 
tories in the perturbed waveguide (See Fig. 
surrounding discussion) . 
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Fig. [17] examines a full-transmission scattering state 
for a large half-stadium with two leads attached at its 
sides. Because scar orbits must self-loop but be other- 
wise unstable[lj, scar states can only participate in trans- 
port when the leads attach at points that are slightly 
displaced from one of the orbit's reflection points; other- 
wise, the classical orbit leaks out the system too quickly. 
The wavefunction in Fig. [IT] shows strong scarring, and 
the multi-modal analysis corroborates the scarring with 
an identifiable classical orbit which just misses the leads. 

Like the square device with obstructions in Fig. [15] 
flux also occurs most strongly along a narrow conduc- 
tance pathway which, in this case, flows along the bot- 
tom of the device while deviating into the bulk at its 
middle. In addition, flux vortices occur throughout the 
system, making interpretation difficult without applying 
our methods. Unlike the square device, however, these 
vortices no longer form identifiable pairs. In the stadium 
state, classical paths do not intersect at 90° angles, but 
take on a variety of other oblique angles. As a result, the 
vortices take on forms that are consistent with the multi- 
modal Husimi map at each intersection. For instance, in 
the black inset, there is strong flow from bottom-left to 



13 



(a) 



•V. 

, iT.m 



Si 





•AC*.*- - 



Figure 17: A scattering wavefunction associated with full 
transmission through the half-stadium in the wavefunction 
representation (a), the Husimi flux (b) and multi- modal anal- 
ysis (c) for a coherent state spread of Ak/k = 10%, indicated 
by the double arrows. The traditional flux from the part of 
the system indicated by the black and red squares is magnified 
in the insets. 



able for informing a design principle in quantum systems, 
since it provides a map of ho w boun daries affect individ- 
ual quantum states (Section |III C ), as well as the im- 
pact of_jx)tentials (Section [4| and magnetic fields (Sec- 



tion IIIIB). Finally, we have shown its utility for illumi- 
nating the many phenomena underlying resonance when 
a closed system interacts with an environment (Section 
IV A), while helping to expl ain th e presence and proper- 
ties of flux vortices (Section [iVB ). Because of its ability 
to contextualize the the flux operator and identify the pri- 
mary conductance pathway in large systems, the Husimi 
projection is an ideal tool for interpreting quantum con- 
ductance simulations. 

This paper focuses on two-dimensional systems, since 
they are ideal for demonstrating the significant physical 
intuition that the Husimi is able to provide. However, its 
definition is not limited to such systems. It is equally 
well suited to three-dimensional systems, and may be 
able to provide a significant contribution to interpret- 
ing molecular orbitals, augmenting such technologies as 
Bader surfaces analysis [26J and local currents [27 \. 



Appendix A: Deriving the Expectation Value of the 

Flux 



In this appendix, we show how to derive Eq. [12] from 
the eigenvalue equation (Eq. [4|. We begin by labeling 
the excited states of the harmonic oscillator at position 
ro oriented along the i th spatial dimension 
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upper-right, with other near- vertical flows forming vor- 
tices, and in the red inset, there are three primary flows 
propagating at 60° to each other, forming the triangular 
arrangement of vortices shown. 



; (ai) 

where e$ is the unit vector along the i th spatial dimen- 
sion. These states form a complete set in which the flux 
operator can be explicitly expressed, using a zero-indexed 
Hermitian matrix, as 



V. CONCLUSIONS 

We have extended the definition of the flux operator us- 
ing a new method for visualizing and analyzing the quan- 
tum wavefunction. By using coherent states as a mea- 
surement operator built off of the Husimi project ion [2 J 
and adapting this technique to generate a vector field, 
we have shown it is identical to th e flu x operator for 
infinitesimal coherent states (Section II B). For finite co- 
herent state spreads, the Husimi projection provides a 
breakdown of the flux into contributions from dominant 
classical paths, which can be retrieved by processing our 
results (Section |IID ). This technique has proven invalu- 
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where A = X a ,i,-\- — There are additional sets of 

harmonic oscillators orthogonal to the above states which 
are centered at points other than r also with zero com- 
ponents in the flux matrix 

The complete set of eigenstates | Ai) , | A2) , | A3) , . . . of 
the flux operator expressed in terms of excited states of 
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the harmonic oscillator are 















(A3) 



with eigenvalues —A, A, and 0. Measurement by the flux 
operator collapses the wavefunction onto one of these 
eigenstates, the infinite majority of which are in the de- 
generate zero-eigenvalue subspace spanning all excited 
states of the harmonic oscillator above |1). Only the first 
two eigenstates, confirmed in Eq. [TT] yield non-zero flux 
values, which, as we have already shown, tend towards 
positive and negative infinity as a — > 0. 

When expanding the flux expectation value, we can 
use the complete eigenbasis to show that 
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From Eq. [TT] it can be shown that the contributions from 
1(^1 0)| 2 and |(^| 1)| 2 cancel themselves due to the op- 
posite sign of the eigenvalues, and only the cross-term 
(^1 0)* (^| 1) — (^| 0) (^1 1)* remains. This form is di- 
rectly related to the commonly-found expression of the 
flux at point rg as 
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Appendix B: Uncertainty Propagation for Husimi 
Vector Addition 



When integrating over the available /c-space in Eq. [24] 
the resulting Husimi flux vector has lower uncertainty 
than the individual terms in the integral, but by how 
much? Understanding this mathematical detail is key 
to appreciating why the Husimi projection is valuable to 
extending the flux operator to an operator with defined 
uncertainty. Moreover, understanding the behavior of 
uncertainty propagation in this integral makes it possible 
to confidently approximate the result with a discrete sum, 
such as the sunbursts in Fig. [2j offering both visual and 
computational advantages. 

We begin by considering the extreme cases. If the 
wavevector orientation remains unchanged for each mea- 
surement, summing up identical measurements has no ef- 
fect on the final relative uncertainty. On the other hand, 
when either the spatial coordinates or the wavevectors 
are sufficiently separated, each Husimi vector constitutes 
an independent measurement; the uncertainty of the re- 
sult will reduce by the square root of the number of mea- 
surements. In general, calculations fall in between these 
two extremes. 



This analysis is concerned with only one dimension, 
since the variance along each orthogonal axis can simply 
be summed. First, the coherent state is expressed in the 
momentum basis as 
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Most generally, the Husimi projection in Eq. [24] is the 
integral of Husimi functions over all of /c-space. In this 
appendix, and in the figures throughout this paper, the 
integral is replaced with a finite sum of test wavevectors 
{k^} which satisfy the dispersion relation at a particular 
energy. 



The variance of the integral in Eq. 24 can be obtained 
by building on intuition about coherent states. It is well- 
known that the /c-space variance of the coherent state 
can be simply derived by integrating the coherent state 
probability amplitude over /c-space, weighting the inte- 
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grand by (k — ko) . Using the notation in Eq. 
gives g\ = yielding the familiar relation a x (Jk 
This can be thought of in the Husimi formulation as a 
statistical result where the quantity is the variance 
of each individual term in the Husimi vector summation. 
In this formulation, the variable is the wavevector and 
the probability function is the probability amplitude of 
the coherent state. Because the probability function is 
complex, we have to take the absolute sum squared. 

Factoring in more than one Husimi function into the 
Husimi projection results in the expression 




2 (k-ki) 2 +i(k-ki)x 



dk, 



(B2) 

where the set {k^} are the set of test wavevectors, pro- 
jected onto the given axis, xq is the spatial point being 
tested, and a is the chosen spatial Gaussian spread. Set- 
ting the coherent states to the same phase at their cen- 
ters, xo — 0, and the above integral can be evaluated to 
return 
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Already it is possible test this result against intuition. 
If each wavevector is identical, then ki — kj =0 and 
the sum of N measurements results in the uncertainty 
which would provide no reduction of relative 
uncertainty. For large values of \ki — kj\ ^> a, the expo- 
nential term will overwhelm the quadratic term and the 
uncertainty becomes o\ = 
uncertainty of y/N. 
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Perhaps most surprising about Eq. |B3| is that the sec- 
ond term, which quantifies the covariance between the 
two measurements, can actually be negative. What are 
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Figure 18: The second term in Eq. |B3| is plotted for the addi- 
tion of two vectors in the Husimi projection. This term repre- 
sents the covariance between the two vectors, and is bounded 
above by 2 and below by —^72 f° r an choices of a. 
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Figure 19: The uncertainty that results from summing three 
vectors of a Husimi projection, as written in Eq. |B3| is plotted. 
The uncertainty is bounded above by 9 /4a 2 and below by 
~ 1.017/4cr 2 . The dashed white lines indicate local minima 
that result from spacing each pair of vectors by \f%j(j, which 
would give a minimum uncertainty for two-vector addition 
(see Fig.[l8). 



its bounds? Fig. 18 plots the quantity Q(k\,k 2 ,(j) = 
2e"^"( /C2-/ci)2 (l - a 2 (k 2 - fci)), showing that a mini- 
mum value of —^372 ~ —0.893 is achieved at \k 2 — k\\ = 
\/3/cr. Every value of \k 2 — k\\ beyond which Q goes 
through zero has achieved nearly independent measure- 
ments, which is found at \k 2 — k\ \ = a -1 . 



The terms in Eq. B3 suggest that when more and more 
vectors are added the uncertainty can be reduced arbi- 
trarily by setting the correct separations between the test 
wavevectors. It even suggests that for three or more vec- 
tors we could possibly produce results with negative un- 
certainty, but intuitively that cannot be possible. To 
appreciate why from an analytical perspective, Fig. [19] 
plots the results of a\ for the addition of three wavevec- 



tors. The minima that occur from maximizing the sep- 
aration between each pair of wavevectors is indicated by 
the white dashed lines. At the center of the graph, a peak 
exists at o\ — 9/4<j 2 , which falls to 3/4<j 2 for areas be- 
yond the area bounded by the white dashed lines, consis- 
tent with earlier observations. There is also a minimum 
(positive) uncertainty which arises from the fact that the 
separation between all pairs of points on a line cannot 



For two vectors the minimum occurs at a 2 
for three o\ w 1.017/4a 2 and for four o\ 



be equal. In Fig. [19] this is evidenced by the fact that 
there are no points where three dashed lines intersect. 

0.981/4a 2 , 

1.017/4cr 2 and for four a\ « 1.036/4a 2 . 
We can generalize and state that for 7V min vectors that 
fall on separate minima, the uncertainty of their sum will 

Moreover, even if vectors are added that do not fall 



on the uncertainty minima in Figs. [18] and [19] they will 
have a negligible impact on the total relative uncertainty. 
So no matter how many vectors contribute to the sum, 
only the vectors on the minima will reduce the relative 
uncertainty, making the key quantity not the total num- 
ber of vectors that are added, but the number that have 
sufficient separation to fall on the uncertainty minima. 

How many vectors is this? We know, for instance, 
that this minimum occurs when the maximum number 
of vector pairs have a separation near and that 

this is likely to occur when they are evenly spaced on a 
line at that separation. Thus we propose that the num- 
ber of vectors that can fall on the minima is given by 
N min = floor (2k(E)a/V3), and using hk = V2mE, we 

can rewrite this as N m i n = floor (cr^J • Substituting 
this value results in the proportionalities 



Ak/k oc 



(B4) 



This makes sense intuitively: the relative uncertainty of 
a finely sampled Husimi vector addition goes down with 
larger a and energy. 

This result deepens the connection between the flux op- 
erator and the Husimi function for small <r, since for very 
small coherent states, the uncertainty minima, which are 
separated by a -1 , grow increasingly far apart. There is 
only a finite range of wavevectors which satisfy the dis- 
persion relation at a given energy, meaning that as the 
coherent states get smaller, fewer and fewer samples in 
fc-space minimize the uncertainty. In fact, at the ex- 
treme limit of a — y 0, the uncertainty cannot be mini- 
mized beyond a single measurement in each orthogonal 
direction, indicating that results for these small coherent 
states have undefined uncertainty, just like the flux op- 
erator. We corroborate this result using a different proof 



in Eq. 20 
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Appendix C: The Hamiltonian 

Numerical simulations in this paper use a free-particle 

2 

Hamiltonian H = — + U(r) sampled on a square 
grid with spacing a and where U(r) = at all points 
unless otherwise stated. This Hamiltonian can be ex- 
pressed in more familiar language by using the tight- 
binding approximation. In this approximation, the effec- 
tive mass envelope function Hamiltonian becomes H = 
e^aja^ — tY2uj\ a ! a j where a^ is the annihilation op- 
erator for the i th grid point, is the energy of the sys- 
tem plus the disorder potential, and the set (ij) cycles 
through all nearest-neighbor pairs. This gives the hop- 



ping term t 



2ma 2 

where zero unless otherwise stated 



and €i = At + Ui, where Ui is every- 



Sec. IIIB 



uses the Peierls substitution [28J to incor- 
porate magnetic fields, using the language of the tight- 
binding model. In this model, the magnetic field con- 
tributes a phase to the hopping potential t: 



tij = t exp [i<j)] , <j) = qA • (r^ — Yj)/ h, 



(CI) 



where is the position vectors of the site corresponding 
to the i th column of the Hamiltonian, h is Planck's con- 
stant, and q is the electron charge. Calculations in this 
paper assume that the magnetic field is perpendicular to 
the plane on which the system sits and is bounded by a 
cylinder centered on the system's center. The radius of 
this column is chosen to be greater than the size of the 
system. Accordingly, the gauge of the magnetic potential 
for an out-of-plane magnetic field is defined such that 



A(r) = — / B z dxdy, 



/■ 



(C2) 



where the integral is over a disc centered on the origin 
and limited by radius r. 

The cyclotron radius can be determined by the relation 



hk 
V ~B^q' 

For a free particle, hk = \/2mE, giving 



This means that at E 
Sec. 



V2mE 
B^qa 2 

: 0.2 



(C3) 



(C4) 



IIIB 



ma 2 



the energy used in 



a magnetic field strength of Bq = 2 x 10 3 -\ 
is sufficient to produce a cyclotron radius that is 2/3 of 
the system radius. This relation is used to predict the 
cyclotron radius for all calculations in this paper. 



conductance calculation. The numerical Green's function 
formalism is used to obtain the scattering wavefunction 
for these calculations, for which modern implementations 
are outlined in several texts|29-31j . In this formalism, 
the Hamiltonian is divided into a left-lead, central region, 
and right-lead projections 



H L Vlc 
H = I V ] LC H c V RC 
tft H R 



(Dl) 



The semi-infinite Green's function at the surface of each 
lead is calculated using the Lopez-Sanchos method, writ- 
ten it as gL,R,(E) fc> r the left (L) and right (R) leads, 
which are both equal. To compute the Green's func- 
tion for the infinite system within the device G(E) the 
semi-infinite surface Green's functions g(E) for each 
lead [32 , 33J are first computed and matched to the surface 
Green's function of the device region, using the numerical 
technique outlined in Mason et al.[34j. 

The coupling matrix for the left lead to the central re- 
gion is then defined by T L (E) = 21m \vl c g L (E)V LC \. 
This results in a density matrix of coherent scattering 
wavefunctions p = GT^G where we have dropped the 
implicit energy dependence. Each coherent scattering 
wavefunction in the system can be obtained by diago- 
nalizing p. Associated with each eigenvector of p will be 
an eigenvalue equal to the likelihood of measuring the 
wavefunction within the system. Since there are gener- 
ally more basis sets within the central region than modes 
available to the system through the semi- infinite leads, 
the vast majority of the eigenvalues will be zero, and 
the number of non-zero eigenvalues will be equal to the 
number of modes available to the system at the given 
energy. This number determines the maximum trans- 
mission across the central region. 

Since a resonant state "traps" the wavefunction at a 
specific energy, it cre ates a striking peak in the den- 
sity of states (see Fig. IV A). As a result, the resonant 
state can be easily identified among the eigenvectors of 
the density matrix since it will be associated with the 
largest eigenvalue near the resonance energy. When dis- 
cussing resonant wavefunctions, it will be assumed that 
we are using a density matrix near the resonance en- 
ergy and examining the eigenvector associated with the 
largest eigenvalue (and measurement probability) at that 
energy. This makes it possible to distinguish the resonant 
wavefunction from other modes which are propagating 
through the system but are unaffected by the resonance. 



Appendix D: Scattering Wavefunctions 

Diagonalizing the Hamiltonian to examine eigenstates 
of a closed system is straightforward. Sec. |IV A[ how- 
ever, examines an open system in a standard ballistic 



Acknowledgments 

This research was conducted with funding from the De- 
partment of Energy Computer Science Graduate Fellow- 
ship program under Contract No. DE-FG02-97ER25308. 



17 



MFB and EJH were supported by the Department of En- ergy, office of basic science (grant DE-FG02-08ER46513). 



[1] Eric J. Heller. Bound-state eigenfunctions of classically 

chaotic hamiltonian systems: Scars of periodic orbits. 

Phys. Rev. Lett, 53(16):1515-1518, Oct 1984. 
[2] K. Husimi. Some formal properties of the density matrix. 

Proc. Phys. Math. Soc. Jpn., 22:264-314, 1940. 
[3] Douglas J. Mason, Mario F. Borunda, and Eric J. Heller. 

A semiclassical interpretation of probability flux. arXiv, 

2012. 

[4] W. Gale, E. Guth, and G. T. Trammell. Determina- 
tion of the quantum state by measurements. Phys. Rev., 
165:1434-1436, Jan 1968. 

[5] Yakir Aharonov and Lev Vaidman. Measurement of the 
schrodinger wave of a single particle. Physics Letters A, 
178:38 - 42, 1993. 

[6] M. Daumer, D. Diirr, S. Goldstein, and N. Zanghi. On 
the quantum probability flux through surfaces. Journal 
of Statistical Physics, 88:967-977, 1997. 

[7] E. J. Heller. Wavepacket dynamics and quantum chaol- 
ogy In M. J. Giannoni, A. Voros, and J. Zinn-Justin, 
editors, Proceedings of the 1989 Les Houches Summer 
School on "Chaos and Quantum Physics", pages 546-663, 
North-Holland, 1989. Elsevier Science Publishers B.V. 

[8] P.W. O'Connor and S. Tomsovic. The unusual nature of 
the quantum baker's transformation. Annals of Physics, 
201(l):218-264, 1991. 

[9] M.S. Child, G. Bruun, and R. Paul. Short time quantum 
phase space dynamics at a 1:2 fermi resonance. Chemical 
Physics, 190:373 - 380, 1995. Overtone Spectroscopy and 
Dynamics. 

[10] Andrea Damascelli. Probing the electronic struc- 
ture of complex systems by arpes. Physica Scripta, 
2004(T109):61, 2004. 

[11] M.C. Gutzwiller. Chaos in classical and quantum me- 
chanics. Interdisciplinary applied mathematics. Springer- 
Verlag, 1990. 

[12] L. Kaplan and E. J. Heller. Measuring scars of periodic 
orbits. Phys. Rev. E, 59(6):6609-6628, Jun 1999. 

[13] W. E. Bies, L. Kaplan, M. R. Haggerty, and E. J. Heller. 
Localization of eigenfunctions in the stadium billiard. 
Phys. Rev. E, 63:066214, May 2001. 

[14] R. W. Robinett. Visualizing the solutions for the cir- 
cular infinite well in quantum and classical mechanics. 
American Journal of Physics, 64(4):440-446, 1996. 

[15] L. A. Bunimovich. On the ergodic properties of some 
billiards. Funct. Anal. App., 8:73-74, 1974. 

[16] M. V. Berry. The bakerian lecture, 1987: Quantum 
chaology. Proceedings of the Royal Society of London. A. 
Mathematical and Physical Sciences, 413(1844) :183-198, 
1987. 

[17] Patrick W. O'Connor and Eric J. Heller. Quantum lo- 
calization for a strongly classically chaotic system. Phys. 
Rev. Lett, 61:2288-2291, Nov 1988. 

[18] Michael V Berry. Quantum chaology, not quantum chaos. 



Physica Scripta, 40(3):335-336, 1989. 

[19] S. Sridhar and E. J. Heller. Physical and numerical ex- 
periments on the wave mechanics of classically chaotic 
systems. Phys. Rev. A, 46:R1728-R1731, Aug 1992. 

[20] Steven Tomsovic and Eric J. Heller. Long-time semiclas- 
sical dynamics of chaos: The stadium billiard. Phys. Rev. 
E, 47:282-299, Jan 1993. 

[21] Fernando P. Simonotti, Eduardo Vergini, and Marcos 
Saraceno. Quantitative study of scars in the boundary 
section of the stadium billiard. Phys. Rev. E, 56:3859- 
3867, Oct 1997. 

[22] Alex Barnett, Doron Cohen, and Eric J. Heller. Deforma- 
tions and dilations of chaotic billiards: Dissipation rate, 
and quasiorthogonality of the boundary wave functions. 
Phys. Rev. Lett, 85:1412-1415, Aug 2000. 

[23] Alex Barnett, Doron Cohen, and Eric J. Heller. Deforma- 
tions and dilations of chaotic billiards: Dissipation rate, 
and quasiorthogonality of the boundary wave functions. 
Phys. Rev. Lett, 85:1412-1415, Aug 2000. 

[24] Doron Cohen, Alex Barnett, and Eric J. Heller. Para- 
metric evolution for a deformed cavity. Phys. Rev. E, 
63:046207, Mar 2001. 

[25] Herman Feshbach. Unified theory of nuclear reactions. 
Annals of Physics, 5(4):357 - 390, 1958. 

[26] R F W Bader. Atoms in Molecules: a Quantum Theory. 
New York: Oxford University Press, 1990. 

[27] Gemma C. Solomon, Carmen Herrmann, Thorsten 
Hansen, Vladimiro Mujica, and Mark A. Ratner. Ex- 
ploring local currents in molecular junctions. Nat Chem, 
2(3):223-228, 03 2010. 

[28] R. Peierls. Zur theorie des diamagnetismus von leitungse- 
lektronen. Zeitschrift fur Physik A Hadrons and Nuclei, 
80:763-791, 1933. 10.1007/BF01342591. 

[29] S. Datta. Electronic Transport in Mesoscopic Systems. 
Cambridge University Press, Cambridge, 1997. 

[30] S. Datta. Quantum transport: atom to transistor. Cam- 
bridge University Press, 2005. 

[31] D.K. Ferry and S.M. Goodnick. Transport in nanos- 
tructures. Cambridge Studies in Semiconductor Physics 
and Microelectronic Engineering. Cambridge University 
Press, 1999. 

[32] M. P. Lopez- Sancho and J. Rubio. Quick iterative scheme 
for the calculation of transfer matrices: application to mo 
(100). J. Phys. F.: Met. Phys., 14:1205-1215, 1984. 

[33] M. P. Lopez-Sancho and J. Rubio. Highly convergent 
schemes for the calculation of bulk and surface green 
functions. J. Phys. F.: Met. Phys., 15:851-858, 1985. 

[34] Douglas J. Mason, David Prendergast, Jeffrey B. Neaton, 
and Eric J. Heller. Algorithm for efficient elastic trans- 
port calculations for arbitrary device geometries. Phys. 
Rev. B, 84:155401, Oct 2011. 



